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ABSTRACT 

Using the integral field spectrograph OSIRIS, on the Keck II telescope, broad near-infrared H and K- 
band spectra of the young exoplanet HR8799b have been obtained. In addition, six new narrow-band 
photometric measurements have been taken across the H and K bands. These data are combined 
with previously published photometry for an analysis of the planet's atmospheric properties. Thick 
photospheric dust cloud opacity is invoked to explain the planet's red near-IR colors and relatively 
smooth near-IR spectrum. Strong water absorption is detected, indicating a Hydrogen-rich atmo- 
sphere. Only weak CH4 absorption is detected at K band, indicating efficient vertical mixing and a 
disequilibrium CO/CH4 ratio at photospheric depths. The i?-band spectrum has a distinct triangular 
shape consistent with low surface gravity. New giant planet atmosphere models are compared to these 
data with best fitting bulk parameters, Tes = HOOK ±100 and \og{g) = 3.5 ± 0.5 (for solar com- 
position). Given the observed luminosity (logLobs/^© ^ ~5.1), these values correspond to a radius 
of 0.75 Rjup 

iai2 mass ^ 0.72 Mjup^Q g - strikingly inconsistent with interior/evolution models. 
Enhanced metallicity (up to 10 x that of the Sun) along with thick clouds and non-equilibrium 
chemistry are likely required to reproduce the complete ensemble of spectroscopic and photometric 
data and the low effective temperatures (< lOOOK) required by the evolution models. 
Subject headings: planetary systems - stars: atmospheres - stars: low-mass, brown dwarfs 



1. INTRODUCTION 

Direct imaging of exoplanets is now becoming a reg- 
ular occurrence with the recent discoveries of 2M1207b 
(Chauvin et al. 2005), Fomalhaut b (Kalas et al. 2008), 
Beta Pic b (Lagrange et al. 2009), IRXS J1609 b 
(Lafreniere et al. 2008), and the quadruple-planetary sys- 
tem HR8799 b, c, d, and e (Marois et al. 2008, 2010). 
HR8799 is, at present, a rare find among exoplanets with 
its four potentially massive planets (~ 5 - 10 Mjup) or- 
biting a ^ 30 Myr A5 star (Zuckerman et al. 2011), and 
is the first and, so far, only multi-planet system to be 
directly imaged. All four planets are accessible to pho- 
tometric and spectroscopic followup observations across 
many wavelengths and, in this paper, near-IR spectra of 
HR8799b are folded into a large set of photometric data 
for a broad-band analysis of the planet's photospheric 
properties. 

Since its discovery, a number of new photometric and 
astrometric measurements of the HR8799 system have 
been made. These new data have broadened the photo- 
metric wavelength coverage and extended the time cov- 
erage of orbital motion (Metchev et al. 2009; Lafreniere 
et al. 2009; Hinz et al. 2010; Currie et al. 2011). The very 
precise astrometric data have already inspired a number 
of papers focused on the orbital dynamics and formation 
of the system; however, the conclusions about the plan- 
etary masses and long term stability of the system from 
these works are mixed (Fabrycky & Murray-Clay 2010; 



Reidemeister et al. 2009; Gozdziewski & Migaszewski 
2009; Dodson- Robinson et al. 2009). As this paper fo- 
cuses on photospheric properties, the dynamical nature 
of the system will not be discussed here. 

The goal of this paper is to characterize the proper- 
ties of HR8799b, independent of what interior/evolution 
models predict and independent of what the age or dy- 
namics of the system suggest. Following a description 
of the new spectra and photometry added by this work 
(Section 2), the available set of photometric and spec- 
troscopic data are compared to similar data for brown 
dwarfs - the closest analogues to young giant planets 
(Section 3). The atmospheric cloud and chemical prop- 
erties of HR8799b are inferred by comparing a new set of 
substellar atmosphere models to the available set of pho- 
tometry and spectroscopy (Sections 4 and 5). The results 
are compared to the predictions of interior/evolution 
models. Also, given the potential for similarities, a com- 
parison is made between HR8799b and 2M1207b (Section 
6). The conclusions from this work are summarized in 
the last section. 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. OSIRIS Spectroscopy 

Observations of HR8799b were made on July 22, 23, 
and 30 (UT) of 2009 and on July 11 and 13 (UT) of 2010 
at the W. M. Keck Observatory. Natural guide star adap- 
tive optics and the OH-Suppressing Infra-Red Imaging 



2 



Barman et al. 



Table 1. OSIRIS Observations 



target 


#exp 


exp. time (s) 


band pass UT date 


1 1 IVO 1 tyif U 


9 


900 


K 


07-22-09 


HD210501 


5 


30 






Rn-l-14 4774 


4 


30 






HR8799b 


6 


900 


H 


07-23-09 






30 






BD+14 4774 


5 


30 






HR8799b 


7 


900 


H 


07-30-09 


BD+14 4774 


7 


30 






HR8799b 


9 


900 


K 


07-11-10 


HD210501 


3 


30 






HR8799b 


8 


900 


H 


07-13-10 


HD210501 


3 


30 







Spectrograph (OSIRIS) instrument (Larkin et al. 2006) 
were used on the Keck II telescope. Broacl-band H and 
K spectra were obtained over a 0.32 x 1.28" patch of 
the sky surrounding the planet at 0.02" spaxel"^ (see 
Fig. 1), where spaxel refers to one OSIRIS spatial reso- 
lution clement. OSIRIS is an integral field spectrograph 
that uses a grid of lenslets in the focal plane to dissect 
the field of view (FOV) . The light from each lenslet is dis- 
persed by a grating, producing an array of spectra that 
are reassembled in the data reduction pipeline (DRP) 
to a [x, y, A] data cube. A scries of 21 iJ-band and 18 
K-hand 900 second exposures were obtained. Appropri- 
ate dark, sky, and AO telluric standard star observations 
were also obtained. See Table 1 for a summary of these 
observations. 

The OSIRIS DRP was used to produce basic calibrated 
data (BCD), and includes spatial rectification of the raw 
data, wavelength and dispersion solutions, sky subtrac- 
tion, cosmic ray rejection, and so forth (Krabbe et al. 
2004). On all three nights in 2009, OSIRIS was operat- 
ing above its normal detector temperature (~ 68K). On 
July 22, the temperature was ~ 79K, but holding steady. 
On July 23, the temperature was ~ 80K and also stable 
throughout the night. However, on July 30, the tempera- 
ture was not stable and increased from 83.56K at the be- 
ginning of the observations to 84.54K at the end. During 
this period of unstable temperatures, the detector noise 
in the dark frames increased exponentially. To avoid tak- 
ing 900 sec sky and dark frames repeatedly throughout 
the night, "scaled darks" were generated for each 900 sec 
exposure on the science target by scaling actual darks 
taken at the beginning and end of the half- night. An 
exponential function was fit to the total counts in each 
actual dark from all three nights and, using this function, 
scaled darks were generated for the science-target cubes 
on July 30. In 2010, OSIRIS was functioning perfectly, 
and the observing conditions were excellent. 

The final BCD data cubes were, at every spatial posi- 
tion, divided by the spectrum of an AO standard star and 
multiplied by a black body (T = 9600K). The spectrum 
of the AO star was extracted from its data cube by fit- 
ting a 2D-Gaussian at every monochromatic image-slice, 
rather than using the fixed-aperture extraction routine 
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Fig. 1. — An unprocessed NIRC2 image (20 second exposure) 
of HR8799. The rectangle shows the orientation and size of the 
OSIRIS spectrograph FOV (0.02" scale). HR8799b is barely visible 
in the center of this rectangle and comparable in brightness to the 
speckle contamination. Dithering was done along the long-axis of 
the FOV. 



in the OSIRIS DRR BD+14 4774 was discovered to be a 
binary (0.07" separation) with a likely mid to late-type 
companion; for this star a two-component Gaussian PSF 
was used to extract the AO fluxes. The final BCD are 
cubes with wavelength and two spatial dimensions. At 
both H and K bands, the planet was clearly visible in 
all median-collapsed cubes (median along wavelength). 
The nominal resolving power of OSIRIS is i? ~ 4000 
with 1600 wavelength channels. However, to improve the 
signal-to-noise ratio, the BCD cubes were binned to 25 
wavelength channels, down to i? ~ 60. The wavelength 
positions of OH sky-lines and telluric star features were 
measured and no artificial wavelength shifts were discov- 
ered; a potential concern given the temperature varia- 
tions encountered during the observations ^. 

While the planet is easily detected in all of the BCD 
cubes, contamination by scattered star light remains a 
concern. Figure 1 shows the FOV and orientation for the 
OSIRIS observations and the extent of speckle contam- 
ination in the observed region. Since the spatial scaling 
of the speckles, radially from the star, is proportional 
to A, speckles can intersect the spatial location of the 
planet in a wavelength-dependent manner. To suppress 
the speckle pattern, each BCD cube was processed with 
custom IDL routines that first subtracts a (spatially) low- 
pass filtered version of the image slice in each wavelength 
channel and robins the data cube (by taking the median 
in the A dimension) to 25 wavelength channels (R ~ 60). 
Each monochromatic image-slice of this rebinned cube is 
magnified about the stars position by A^/A, where A^ is 
the median wavelength in either H or K band. In these 
shifted data cubes, individual speckles align and arc fit 
with simple polynomial functions of A. Prior to this step, 

^ Several updates of the OSIRIS DRP became avail- 
able during the course of this work; one of these was 
directly related to preventing unwanted wavelength shifts 
(http : // irlab . astro . ucla . edu/osir is/). 
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Fig. 2. — Left: A single H-band basic calibrated data cube 
(900s exposure at one dither position), median combined along 
the wavelength dimension. The planet is visible in the upper left 
quadrant. Bright speckles are visible across most of the image. The 
image is oriented such that the star is to the right. Right: Same 
data as in the left panel, but after speckle-suppression. 



the planet is masked to avoid biasing the polynomial fit. 
The final polynomial fits to the speckle fluxes are sub- 
tracted from the shifted cube before finally reversing the 
magnification of the cube, placing the planet back in its 
original position. Collapsed i?-band data cubes, before 
and after speckle removal, are shown in Fig. 2 (a simi- 
lar improvement is seen in the before and after ii'-band 
data cubes). The speckle suppression scheme described 
here is similar to that used for the low-mass companion 
GQ Lupi b (McElwain et al. 2007); however, McElwain 
et al. did not rescale the cubes for speckle fitting but 
instead fit the underlying speckle halo in each spectral 
bin of the data cube since each spatial location is more 
than critically sampled. 

The speckle suppression steps described above contain 
two important variables: the size of the masked region 
(centered on the planet) and the order of the polynomial 
used to fit the speckle fluxes. To determine the best pair 
of parameters, ten fake planets (similar in brightness and 
FWHM as the real planet but with constant F^) were 
inserted into the BCD data cubes, evenly spaced but 
avoiding the real planet location as well as the edges of 
the FOV. These cubes were processed using the twelve 
combinations of four mask radius values (from 0.5 to 
2 in 0.5 spaxel steps) and three polynomial orders (1, 
2 and 3). The parameters resulting in extracted fake 
planet spectra with the smallest mean deviation from 
a flat line and smallest RMS(A) were mask radii equal 
to 1.5 and 2.0 for H and K bands respectively and a 
polynomial order equal to 1 for both bands. The mean 
residuals of all ten fake planet spectra are shown in Fig. 
3, straddling zero flux. The RMS (A) for the fake planet 
spectra were comparable to, or smaller than, the final 
uncertainties adopted for the planet spectrum, provid- 
ing confidence that the errors have not been significantly 
underestimated. Furthermore, the data collected in 2009 
compare very well to those from 2010, indicating that the 



thermal instabilities and scaled darks had little impact 
on the final reduced spectra. 

The planet fluxes were obtained from the speckle- 
suppressed data (SSD) cube by convolving each 
monochromatic image-slice with a two-dimensional, 
Gaussian, point-spread-function with FWHM = 2.2 {H- 
band) and 2.7 (-fC-band) spaxels. The final spectra were 
flux calibrated using i?-band (Metchev et al. 2009) and a 
revised Ks magnitudes (discussed below). The final ex- 
tracted spectra and corresponding error-bars (taken as 
the RMS of all the extracted spectra in each band-pass) 
are shown in Fig. 3 and listed in Table 2. 

As can be seen from the residuals of the fake planet 
fluxes (Fig. 3), the noise in the spectra is correlated 
between different wavelengths. This is a natural con- 
sequence of the fact that the noise is primarily from the 
speckle pattern - the typical speckle has a monochromatic 
size of A/D, as does the extraction box for the plane- 
tary spectrum. A single speckle's position scales with 
wavelength proportional to A. At K band, the planet 
is located approximately 42 A/D from the star. For the 
speckle to move completely across the extraction box re- 
quires a change in wavelength of (43/42-1) x 2.1 fim = 
0.05 fiui. Hence speckle noise patterns will be correlated 
across the spectrum over a range of 0.05 /im (Sparks 
& Ford 2002). In principal this should be formally ac- 
counted for in the estimate of when fitting models to 
the data (described below). However, since the spectral 
features being fit are also of width ^0.05 or larger, 
the correlations in the noise have been ignored. 

The planet spectra, before and after speckle suppres- 
sion, are also compared in Figure 3. Near 1.7 /im, the 
impact of a bright speckle can be seen as a fairly sharp 
peak spanning two wavelength channels. Speckle noise 
and photon noise are roughly equal in these data sets at 
this binning, so speckle suppression improved the SNR 
by approximately 50%. Most importantly, the speckle 
noise that was removed was spectrally correlated across 
multiple channels and non-Gaussian, e.g. the feature at 
1.7 /im in Figure 3, and would have a much more sig- 
nificant effect on the final interpretation than the white 
photon noise. 

2.2. Photometry 

On August 7 2009 (UT), when OSIRIS was too warm 
to operate, an angular differential imaging (ADI, Marois 
et al. 2006) sequence was taken with NIRC2 using 
four narrow-band filters in the K-hand window (see 
Table ??). These filters were selected to sample, as 
uniformly as possible, the full K-hand region accessi- 
ble with OSIRIS and thus serve as a consistency check 
on the shape of the final, fully processed, OSIRIS 
spectrum. These narrow-band data are compared in 
Figure 3 to equivalent flux points obtained from the 
OSIRIS data convolved with the filter response profiles. 
The NIRC2 narrow-band colors are consistent with the 
OSIRIS spectral shape on either side of the peak fiux; 
however, HelB (2.06 /im) photometry is only consistent 
at ^ 2-CT. The overall agreement, between these two 
non-contemporaneous data sets, processed with differ- 
ent reduction software, provides some confidence that 
the planet's spectrum was faithfully extracted from the 
speckle-contaminated data. 

Also in 2009, new CH4 (4% S/L) photometry was 
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Fig. 3. — OSIRIS H and K-hand fluxes of HR8799b (scaled to 10 parsecs) plotted with 1-cr uncertainties. The location of prominent 
water, CH4, and CO absorption bands are indicated. The fluxes extracted without the speckle suppression algorithm are shown as dotted 
lines. The bottom two curves are the mean residuals of fake planets with flat spectra extracted from the same data cubes (see text for 
details). The Kn3 spectrum of Bowler et al. (2010) is shown as green pluses (scaled arbitrarily down) over plotted with the broad-band 
spectrum (black points) interpolated onto the Bowler et al. Kn3 wavelength points. The mean l-cr uncertainties across the Kn3 range are 
shown at either end for each data set. The larger red filled symbols are the NICI CH4 short/long (boxes), NIRC2 narrow (circles), and 
NIRC2 broad-band (stars) photometry taken from Table ??. Blue symbols are the corresponding photometry derived from the OSIRIS 
spectra. 



Table 2. OSIRIS H and K Spectra (scaled to 10 parsecs) 



A Fi, (mjy) error (1-cr) A Fi, (mjy) error (1-cr) 



obtained with the Near Infrared Coronagraphic Imager 
(NICI; Toomey et al. 2003 ) on Gemini South. The NICI 
CII4 short/long filters more optimally probe the strength 
of CH4 absorption than the NIRC2 CH4 short /long fil- 
ters. The former are narrower, have less overlap, and 
have central wavelengths that correspond well to the min 
and max fluxes seen in H-hand spectra of mid to late 
T dwarfs. A comparison between these new photomet- 
ric data and equivalent band-integrated flux points from 
the OSIRIS spectrum is shown in Fig. 3. Here again the 
CII4 short/long slope is in excellent agreement with the 
OSIRIS spectrum. 

Since the discovery of the IIR8799 planetary system 
more broad-band observations have been published and 
improvements have been made in the LOCI (Lafreniere 
et al. 2007) and ADI algorithms (at least those used by 
this group). For the analysis presented here, the i7-band 
photometry of Metchev et al. (2009) has been adopted. 
Also, Ks data taken in 2010 have been analyzed and a 
new Ks magnitude (Mks — 14.15 ±0.1) has been ob- 
tained. See Table ?? for a full list of the NIRC2 and 
NICI photometry used in this study. 

2.3. Comparison to KnS spectrum 

On June 21st (UT) 2009, Bowler et al. (2010) observed 
HR8799b using the OSIRIS Kn3 narrow-band filter. A 
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Table 3. NIRC2 and NICI Photometry 



Filter Ac Mag Fiopc (mjy) ZP(Jy) Rcf 



NIRC2 broad-band 

1 05 16.30 0.46 ± 0.07 1521.1 Mar08 

H 1.63 15.08 0.94 ± 0.11 1010.0 Met09 

Ka 2.15 14.15 1.43 ± 0.11 654.2 updated 

L' 3.78 12.66 2.06 ± 0.21 238.5 Mar08 

NICI narrow-band 

Ca^^short 135 15.18 0.88 ± 0.14 1035.7 new 

CH4,o„g 1.68 14.89 1.06 ± 0.18 959.7 new 

NIRC2 narrow-band 

He 1 H Tm 14.73 0.85 ± 0.07 665.9 new 

H2{v=l-0) 2.13 14.16 1.52 ± 0.11 703.4 new 

BrT,2 2.16 13.93 1.66 ± 0.12 618.9 new 

H2(v=2-1) 2.26 14.09 1.37 ± 0.11 593.5 new 




1.6 l.S 2.0 2.2 2.4 

Wavelength (^m) 

Fig. 4.— OSIRIS H anad X-band spectra of HR8799b (plotted 
as Fx) compared to field brown dwarf spectra (taken from the 
SpeX spectral library). Each band was normalized individually; 
the actual J-Ks color is given below each object name. All of 
these field objects are bluer than HR8799b. 



over 300 low-resolution L and T dwarf spectra from ~ 
0.65 to 2.55 /zm and includes a number of very red T 
dwarfs. These spectra were rebinned to match the reso- 
lution and sampling of the HR8799b OSIRIS spectra and 
used to find the best match to the H and K band spec- 
tra separately by minimizing x^- All L and T dwarfs 
available in the SpeX library were used, except known 
binaries. 

Figure 5 shows the results of fitting the SpeX observa- 
tions to HR8799b. There is a clear minima between 
TO and T2, with the best fitting objects having spectral 
types T1.5 and TO. 5 for the H and K bands, respectively 
(within the range found by Bowler et al., who fit their 
Kn3 spectrum to the same SpcX library). An empiri- 
cal fit to the broad-band photometry was also performed 
by Bowler et al.; their best overall match (a peculiar L6 
2MASS 2148-f4003) is compared to the H and /C-band 
spectra in Fig. 5. While this L6 is a good match to the 
broad band photometry, it is not a good match to the 
broad-band near-IR spectral data. 

It is interesting to note that, while fitting brown dwarf 
spectra is completely independent of how well the broad- 
band colors agree (each band was normalized individu- 
ally for the calculation), there is a clear correlation 
between J-Kg color and x^. At first glance, this is simply 
restating that the best fits are found among the late L 
and early T dwarfs, which have the reddest J-Ks colors. 
However, in this case, the best fitting objects are notice- 
ably redder than the mean J-Kg color for their spectral 
type. The lower panel of Fig. 5 shows versus J-Ks 
color for all T dwarfs. For both bands, the best fitting 
objects (a TO. 5 and a Tl) are redder than the mean J- 
Ks for TO to T2 (vertical dotted lines, values taken from 



Mar08 = Marois et al. (2008); Met09 = Metchev et al (2009) 



comparison between the narrow and broad-band spectra 
is shown in Fig. 3, and both spectra agree within the 
l-CT error-bars. The Bolwer et al. spectrum decreases 
in fiux more noticeably than the broad-band spectrum; 
Bolwer et al. attribute this decrease in fiux to possible 
weak CII4 absorption. The broad-band spectrum, which 
encompasses much more of the CH4 absorption band, 
indicates an even weaker CH4 signature than indicated 
by the Kn3 spectrum. 

3. EMPIRICAL COMPARISONS 

The H and if-band spectra of IIR8799b show several 
interesting features. In the H band, a pronounced trian- 
gular shape indicative of weak collision induced absorp- 
tion (CIA) and low surface gravity is seen. The spec- 
trum also shows no evidence of strong methane absorp- 
tion (consistent with CII4 on/off photometry) as would 
otherwise be expected in cold T-type brown dwarfs. At 
K band, the spectrum shows very deep water absorption 
bands. As in the H band, there is no evidence of strong 
methane absorption or very strong CO absorption. 

HR8799b occupies sparsely populated regions of near- 
IR color-magnitude diagrams. Despite this fact, it is still 
useful to compare the observed spectra of IIR8799b to 
those of known brown dwarfs. Figure 4 compares the H 
and iiT-band spectra to a sequence of L and T dwarf spec- 
tra from type T6 to LI. The near-IR spectrum HR8799b 
is fairly distinct from these typical field brown dwarfs, 
which generally have weaker water absorption and, for 
the T dwarfs, have deeper CII4 absorption than seen 
in HR8799b. The fact that none of the hotter, L-type, 
brown dwarfs provide a good match to HR8799b is con- 
sistent with the low effective temperature deduced from 
cooling track models (Marois et al. 2008). 

Brown dwarf near-IR colors exhibit considerable 
spread within a given spectral type and a number of 
peculiar (anomalously red or blue) L and T dwarfs are 
turning up (Kirkpatrick et al. 2008; Burningham et al. 
2010) . Consequently, while standard BDs do not match 
IIR8799b very well, perhaps some of the more peculiar 
BDs might. The SpeX Prism spectral library^ provides 

^ http:/ /web. mit.edu/ajb/www/browndwarfs/spexprism 
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Fig. 5. — top panel: real brown dwarf spectra that best fit the H 
and K bands of HR8799b separately (red: best fit to H band, blue: 
best fit to K band). Each spectrum is scaled arbitrarily to have 
the same total normalized fiux. Plotted in green is the spectrum 
of 2MASS J2148+4003, the object found by Bowler et al. (2010) 
to be the best overall match to HR8799b. middle panel: values 
versus spectral type for real brown dwarfs compared to H (blue) 
and K (red) bands, bottom panel: values versus J — Kg color for 
real T-type brown dwarfs compared to the H (blue) and K (red) 
bands. Vertical dotted lines indicate the J — Ks color for HR8799b 
as well as typical colors for field TO, Tl, and T2 dwarfs. 



Faherty et al. 2009). While none of the objects available 
in the SpeX library are as red as HR8799b, there is a 
clear indication that early T dwarfs, with anomalously 
red near-IR colors, are the best match. This trend was 
also independently noted by Bowler et al.. 

The best fit for H band, SDSS J151643.014-305344.4, 
(here after SDSS 1516-^30), has J-K, = 1.77 and has 
been classified as TO. 5 ±1 (Chiu et al. 2006; Burgasser 
et al. 2006) and T1.5 ±2 (Burgasser et al. 2010). While 
SDSS 1516-1-30 matches HR8799b in H band very well, 
the comparison is poor at K band. Stephens et al. (2009), 
fitting their model atmospheres to available near-IR and 
mid-IR spectra of SDSS 1516-^30, found Tes = 1000 to 
1100 K, log(5) = 4.5 and a fairly cloudy atmosphere 
for such a low Tog (fsed = 1 to 2, in the Ackerman & 
Marley (2001) parlance). This Tcfi range is lower than 
typical for a T0.5, by 200 to 300K (Stephens et al. 2009; 
Golimowski et al. 2004). Burgasser et al. (2010) identify 
SDSS 1516-1-30 as a weak, if not unlikely, binary candi- 
date; however, binarity has yet to be tested with high 
resolution imaging. Leggett et al. (2007) have also noted 
this object's mid-IR colors are best reproduced by mod- 
els with thick clouds and stronger CO absorption than 
predicted by pure chemical equilibrium, consistent with 
the high eddy diffusion coefficient used by Stephens et 
al. in their model comparison. 

The best fit for K band is 2MASS J21392676-I-0220226 
(here after 2M2139). The brown dwarf has J-K^ = 
1.68 and has been classified as TO (Reid et al. 2008), 
T1.5 (Burgasser et al. 2006), and recently as T2.5 ±1 
(Burgasser et al. 2010). Not only does 2M2139 match 
IIR8799b reasonably well in the K band, the comparison 
is also about as good as SDSS 1516-1-30 in the H band. 



Burgasser et al. (2010) suggest that this object is a strong 
binary candidate, though no high-resolution imaging has 
been reported. It is quite possible that the object is 
single and simply non-standard with atmospheric prop- 
erties that are a mix of those found in the cloudy and 
cloud-free brown dwarf photospheres. The 7^-band spec- 
trum of 2M2139 is similar to that of SDSS J0758-F42 (a 
T2), the brown dwarf found by Bowler et al. that best 
matched their Kn3 spectrum. The _ff-band spectrum of 
SDSS J0758-I-42, however, is noticeably bluer than the 
HR8799b OSIRIS iJ-band spectrum. 

In summary, the empirical comparisons suggest a spec- 
tral type of Tl ±1, consistent with the T2 match that 
Bowler et al. (2010) found when comparing their Kn3 
spectrum to the SpeX library. However, it appears as 
though only the reddest T dwarfs in this spectral type 
range provide reasonable matches to the HR8799b spec- 
tra. While a Tl is a close match, the near-IR colors 
of IIR8799b are still well out of the normal range for 
this spectral type (being more consistent with the mid L 
dwarfs; consequently, there are probably significant lim- 
itations to what any near-IR spectral type can tell us 
about the basic properties of IIR8799b. 

4. MODEL ATMOSPHERES 

Recognizing that HR8799b falls outside of the color- 
magnitude range of most known substellar objects, com- 
paring the observed data to synthetic spectra from model 
atmospheres is the next logical step. The PHOENIX model 
atmosphere code (Hauschildt et al. 1997) has been used 
to explore a variety of physical mechanisms that might 
explain the unusually red near-IR colors, low luminos- 
ity, and lack of strong methane bands. The version 
of PHOENIX (vl6) used here has been substantially up- 
dated since the introduction of the PHOENIX "cond" and 
"dusty" brown dwarf atmosphere models (Allard et al. 
2001). Some, but not all, of the updates include the 
replacement of the original Allard chemical equilibrium 
solver with a new, more robust and chemically com- 
plete, equation of state'^, a revised treatment of pure dust 
grains, improved alkali line profiles (Allard et al. 2007), 
and new molecular line lists, most notably for methane 
(Warmbier et al. 2009; Hauschildt et al. 2009) and wa- 
ter (Barber et al. 2006). Updates to the cloud modeling 
and the addition of mixing-induced disequilibrium chem- 
istry are the most relevant here and are described in the 
following sections. 

4.1. A Simple Cloud Model 

Atmosphere models born from the study of brown 
dwarfs are commonly used to estimate the properties of 
giant planets and to estimate detection yields for direct 
imaging surveys, so it is natural to use such models as 
starting points for an analysis of HR8799b. Arguably 
the most complex physical process to include in substel- 
lar atmosphere models is the presence of clouds. The 
term cloud is used here to refer to solid or liquid particles 
suspended in the atmosphere and, depending on the tem- 
peratures and pressures, can include a complex mixture 
of species ranging from iron particles to ices. The major- 
ity of substellar atmosphere models are one-dimensional, 

^ The new equation of state routine is known as the Astrophysi- 
cal Chemical Equilibrium Solver, or ACES (Barman &; Hauschildt, 
in prep). 
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Fig. 6. — Near-IR color-magnitude diagrams for brown dwarfs (Leggett et al. 2002; Knapp ot al. 2004; Casewell et al. 2007) and the 
planets 2M1207b and HR8799bcd. Dotted lines show color-magnitude tracks for chemical equilibrium models of constant T^fj {900K to 
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assuming either plane-parahel or spherical geometry, and 
consequently when clouds are included their spatial prop- 
erties can only be modeled radially. Nevertheless, the 
"patchiness" of clouds can be explored in an ad hoc fash- 
ion via a combination of cloudy and cloud-free models 
(Burgasser et al. 2002) or, more recently and selfcon- 
sistently, within a single one-dimensional model (Marley 
et al. 2010). For an extensive review and comparison 
of theoretical cloud modeling in substellar atmospheric 
conditions, see Helling et al. (2008a). 

The observed luminosity and colors of HR8799b have 
already been shown by Marois et al. (2008) to be incon- 
sistent with those of field brown dwarfs. Therefore, by 
extension, these data are inconsistent with the predic- 
tions of most published substellar atmosphere models, 
as they are often tuned to match existing observations. 
A large number of early brown dwarf models reproduced 
the basic trends of L and T dwarfs by approximating the 
cloudy and cloud-free conditions in a very simple man- 
ner; by simply turning on (cloudy) or off (cloud-free) the 
opacity of dust in the atmosphere at its chemical equilib- 
rium temperature-pressure location (Allard et al. 2001). 
This approach, taken to be the extreme limits of cloud 
influence, is not expected to reproduce every individual 
object well, especially those in the transition region from 
spectral types L to T. In fact, most objects should fall 
between these two limiting cases. In order to match the 
photometric and spectroscopic observations of IIR8799b 
a different approach is required - a cloud model that is 
intermediate to these two limits is needed. 

Recent improvements in the atmospheric modeling of 
transition brown dwarfs have been made, with most ef- 
forts attacking the cloud problem by adding one or more 
new parameters that regulate the cloud thickness, par- 
ticle density, grain size distribution and so forth (Ack- 
erman & Marley 2001; Marley et al. 2002; Tsuji 2005; 



Cooper et al. 2003; Burrows et al. 2006). However, 
see Helling et al. (2008b) for a completely different ap- 
proach. The purpose of this paper is not to develop a 
new complex cloud model, but rather to identify the ma- 
jor physical processes that shape the spectral properties 
of HR8799b. To that end, a simple parameterized cloud 
model is adopted for the analysis presented here. The 
lower boundary (cloud base) is established by the inter- 
section of the temperature-pressure (T-P) profile and the 
relevant chemical equilibrium condensation curve (Pc, 
using pressure as a proxy for height) . The behavior of the 
cloud above Pc is determined by a single free parameter 
(Prnin)- The equilibrium dust concentration is multiplied 
by a function that is 1 for Pgas ^ Pmin and decays ex- 
ponentially for Pgas < Pmin- If Pmin > Pc, then the 
maximum dusty-to-gas ratio is also lowered relative to 
an equilibrium cloud. Examples of vertical cloud prop- 
erties are discussed below in more detail. The particle 
sizes follow a log-normal distribution with a prescribed 
modal size (ao) from 1 to 100 fim that is independent of 
height. This cloud model shares some similarities with 
other parameterized cloud models (Tsuji 2005; Burrows 
et al. 2006); however, rather than accounting for a small 
representative set of grains, all thermodynamically al- 
lowed grains with opacities in the PHOENIX databases are 
included in the total cloud opacity; see Ferguson et al. 
(2005) for a list of grains included in the model. 

Figure 6 shows the location of HR8799b relative to L 
and T brown dwarfs in three near-IR color-magnitude 
diagrams. Synthetic color-magnitude sequences are also 
shown for chemical equilibrium models (K^^ — 0) with 
radius set to that of Jupiter. Each sequence corresponds 
to a different Tefi (with \og{g) = 4.0) and different cloud 
thickness varying from no clouds (at the blue end) to 
maximum cloud thickness (at the red end). These se- 
quences of models are not associated with any evolution- 
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ary cooling tracks since the radius has been arbitrarily 
set to that of Jupiter. Such sequences illustrate the color 
trends produced by changing vertical cloud thickness 
(note this is different from the uniform sinking of clouds). 
Figure 6 shows that the slope in a CMD across the L- 
to-T transition region, including the J-band brightening, 
is reproduced very well by atmospheres that have clouds 
smoothly decreasing in vertical thickness. However, as 
pointed out by Burrows et al. (2006), reproducing color- 
trends does not guarantee good reproduction of near-IR 
spectra; the models shown in Fig. 6 have not been com- 
pared to spectra of field brown dwarfs (an exercise left 
to another paper). 

Figure 6 also shows that very red objects with rel- 
atively low luminosities, like HR8799b, can be repro- 
duced by allowing low Toff models (typically assumed 
to be in the regime of cloud-free photospheres) to have 
clouds extending across their photosphere. The location 
of HR8799b, in all three CMDs is intersected by a model 
with Tcff = lOOOK, \og{g) = 4.0, and a cloud that starts 
to fall off near 1 bar of pressure (though the precise Pmin 
matching the position of HR8799b is slightly different in 
each CMD). Figure 6 also suggests that the cloud prop- 
erties of the HR8799 planets and 2ml 207b are similar 
to the those of late-type L dwarfs. The comparison in 
Fig. 6 does not directly identify a best fit to the ensem- 
ble of data, but simply illustrates the impact that cloud 
thickness has on the near-IR colors of an object with typ- 
ical substellar temperatures and gravities and points to 
a reasonable place to start looking for a good match to 
HR8799b, namely Teft~ lOOOK and log(ff) ~ 4.0. 

4.2. Disequilibrium Chemistry and Vertical Mixing 

The majority of atmosphere models for brown dwarfs 
and exoplanets assume local chemical equilibrium; in 
other words, the mole fractions for each chemical com- 
pound is determined by the local temperature and pres- 
sure, and overall element abundances. For most atmo- 
spheric compounds, this is an excellent approximation 
as the chemical reactions that establish the equilibrium 
mole fractions are rapid compared to the mixing time 
scales (even in the convection zone). However, this is not 
always the case and it has been shown that, for example, 
the reactions that ultimately replace CO with CH4 can 
be very slow in Jovian and brown dwarf photospheres. 
The same can be true for N2 and NH3 . Prinn & Barshay 
(1977) showed that vertical mixing of CO can be impor- 
tant in Jupiter and, later, Fegley & Lodders (1996) pre- 
dicted that this would also be the case for brown dwarfs. 
The spectroscopic detections of CO in several cool T-type 
brown dwarfs (NoU et al. 1997; Saumon et al. 2000, 2006; 
Geballe et al. 2009) along with weaker than predicted 
fluxes from 4.5 fim ground and space-based photom- 
etry (Patten et al. 2006; Leggett et al. 2007) have es- 
tablished disequilibrium CO mole fractions as a common 
occurrence in brown dwarfs. Non-equilibrium chemistry 
is, therefore, a plausible mechanism for the apparently 
weak CH4 absorption in HR8799b. 

Disequilibrium CO/CH4 and N2/NH3 chemistry has 
been added to PHOENIX following Smith (1998). The 
reaction time scales (in seconds) for N2 — ?> NH^ and 
CO CHi are from Lodders & Fegley (2002) and Yung 
et al. (1988), respectively: 
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Fig. 7. — CO (dashed lines) and CH4 (dotted lines) mole frac- 
tions for a cloudy, solar metallicity, model atmospheres with T^fj = 
800K and four different surface gravities (corresponding colors in- 
dicated in the figure legend). Chemical equilibrium mole fractions 
are plotted in the left panel. Non-equilibrium mole fractions are 
plotted in the right panel, where the mole fractions are quenched 
at different depths determined by the intersection of the CO chem- 
ical reaction and mixing time scales (tqo and r^ix > dotted lines) . 
The latter time scale was arbitrarily set to 1 yr. In both panels, 
the vertical dotted line marks the approximate location of the pho- 
tosphere (taken here to be where the local temperature equals the 
effective temperature). 

KN2N(H2) 

o /-81515\ 
KN, = 8.45 X IQ-^'exp i ^ j (2) 

^ NjCO) 

KjN{H)N{H2C0) ^ ' 

where T is the gas temperature in Kelvin, n are the rate 
coefficients in cm? / sec (with K7 interpolated from tabu- 
lated values in Yung et al. 1988) and N(species) is the 
number density of the specified species in cm~^ . 

The mixing time scale {Tmix) is computed using mix- 
ing length theory in the convection zone, and set equal to 
L^g/Kzz in the radiative zone, where the effective length 
scale (Lcff) is determined following the recipe in Smith 
(1998). For the range of models explored here, LcS was 
found to be ~ 0.1 to 0.2Hp, with Hp equal to the pres- 
sure scale height. This range of Leff is comparable to that 
found for Jupiter (Smith 1998; Visscher et al. 2010). The 
major unknown when modeling vertical mixing is the 
eddy diffusion coefficient (K^z) which, as is frequently 
done, is assumed here to be an adjustable parameter typ- 
ically ranging from 100 to 10^ cm^ sec~^. With these 
chemical and mixing time scales, the mole fractions of 
CO/CH4 and N2/NH3 are "quenched" at a specific depth 
in the atmosphere where Tmix = Tchem- Below this point 
in the atmosphere, where Tmix > Tchem, chemical equi- 
librium prevails. This simple quenching model has been 
used in several recent studies of brown dwarfs and giant 
planets (Hubeny & Burrows 2007; Fortney et al. 2008). 

The mixing and chemical time scales alone do not de- 
termine the final quenched mole fractions, the atmo- 
sphere temperatures and pressures still establish these 



The Atmosphere of HR8799 b 



9 



values at the quenching depth. Consequently, effective 
temperature, surface gravity, and mctallicity play an im- 
portant role. Of particular importance is the relative lo- 
cation in the atmosphere of equal CO/CH4 equilibrium 
mole fractions {Pco=CHi) and where Tmix = Tchem {Pq)- 
If Pq > Pco=CHi , the CO/ CH4 ratio can be greater than 
1 in cool substellar atmospheres where CO would other- 
wise be predicted to be orders of magnitude less than 
CH4. 

As shown in Figure 7, the magnitude of the disequilib- 
rium effect is sensitive to surface gravity, as this quantity 
acts to shift the locations of Pq and Pco=CHi ■ As gravity 
decreases, the crossing point of the radiative and mixing 
time scales moves to deeper layers (higher T), while the 
crossing of the equilibrium mole fractions of CO and CH4 
moves outward. For K^^^ values typically used in brown 
dwarf models and at low surface gravities expected for 
young exoplancts. the CO/CH4 ratio can become signif- 
icantly greater than 1 indicating a role-reversal in the 
dominant carbon-bearing species at photospheric depths. 
Such a role reversal has not been explored thoroughly 
in brown dwarfs as their surface gravities and effective 
temperatures almost always place the time scale cross- 
ing point at a pressure lower than that of the equilibrium 
CO/CH4 crossing point, resulting in a CO/CH4 ratio al- 
ways less than 1 (though still greater than equilibrium 
chemistry predictions). Hubeny & Burrows (2007) also 
noted a gravity dependence for non-equilibrium chem- 
istry, but did not explore gravities below \og{g)= 4.5 
nor encountered a CO/CH4 reversal. Since a low sur- 
face gravity and effective temperature are possible out- 
comes for HR8799b, a lack of (or very weak) CH4 could 
be caused by disequilibrium chemistry. 

4.3. Photochemistry 

HR8799b is located a mere 68 AU from its host star 

and receives far too little stellar flux for irradiation to 
have a significant influence on its atmospheric thermal 
structure (at near-IR photospheric depths). However, 
HR8799A is much more luminous (~ 5 L©) and hot- 
ter (TefF= 7500K) than the Sun and the UV flux from 
HR8799A may be sufficient to produce photochemical 
reactions in the atmosphere of HR8799b and certainly 
for the inner three planets. Based on a synthetic spec- 
trum from a stellar model atmosphere calculation tai- 
lored for HR8799A (also using PHOENIX, with atoms and 
ions treated in non-local thermodynamic equilibrium), 
HR8799b is estimated to receive ~ 8 times the fiux that 
Jupiter receives from the Sun for A < 2000A, and about 
10 times as much flux as Neptune for A < 3000A. Given 
the larger atmospheric temperatures and possibly ex- 
treme non-equilibrium CO/CH4 mole fractions compared 
to the giant planets in the Solar System, the photochemi- 
cal reactions in HR8799b will also be quite different from 
those found in the Solar System. Consequently, photo- 
chemistry could impact the model comparisons discussed 
below. Unfortunately, photochemical reactions are not 
included in the version of PHOENIX used here. Future ob- 
servations and modeling will be needed to address this 
possibility. 

5. MODEL COMPARISONS 

The complete set of HR8799b broad-band photome- 
try covers a broader wavelength range than the OSIRIS 
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Fig. 8. — Mean contours obtained by comparing model spec- 
tra to the broad-band photometry (Table 2) and OSIRIS H and 
K spectra (Table 3). The best fitting model has Tcfi= HOOK 
and \og(g)= 3.5. The right ordinate indicates the radii required to 
match the observed luminosity of HR8799b given the effective tem- 
peratures on the left ordinate. Using these radii and the surface 
gravities along the abscissa, lines of constant mass can be plotted 
(dashed lines for 0.3 and 1 Mjup, solid line for 20 Mjup). These 
lines of constant mass axe independent of atmosphere or evolution 
model predictions. Points along the dotted line are consistent with 
a Kelvin-Helmholtz cooling time of 1 Myr. 



spectroscopic observations and, therefore, are likely to 
provide the greatest leverage when estimating effective 
temperature and, potentially when estimating the cloud 
properties. Marois et al. (2008) did not perform a model 
fit to the discovery photometry but, rather, focused on 
the parameters inferred from theoretical cooling tracks 
and the observed luminosity and age range. With new 
and improved photometry available, new model fits to 
these data are now warranted (see Bowler et al. (2010) 
and Currie et al. (2011) for independent model compar- 
isons to similar photometric data and different model at- 
mospheres). Additionally, the near-IR spectra are likely 
to contain the greatest information on surface gravity, 
chemical composition (via broad-band absorption fea- 
tures), and also on the cloud properties. The goal in 
this section is to identify a single model atmosphere that 
best fits all three observational data sets and to under- 
stand what physical mechanisms play dominant roles in 
shaping the overall spectral energy distribution. 

5.1. Solar Abundance Model Comparison 

The cloud and disequilibrium chemistry models de- 
scribed above were incorporated into PHOENIX. A modest 
grid of synthetic spectra was computed with five key pa- 
rameters considered: T^s [= 800K - 1500K, lOOK incre- 
ments], log(g) [= 2.5 5.0. in 0.5 increments and CGS 
units], K^^ = lO'' cm2 s^^, [= 10\ 10^ 2xl0^ 

4x10^ 6xl0^ 8xl0^ 1x10^ dynes cm-^j, a,, [ = 1, 5, 
10, 100 /Ltm]. Solar abundances were adopted for this 
grid (Asplund et al. 2005), but the question of non-solar 
metallicity will be discussed below. 

Given the large number of free parameters and the 
sparse sampling of the model grid, obtaining a traditional 



10 



Barman et al. 





1 100K, 109(g) = 3,5 













2 3 
Wovelength (m^) 




Wovelength (m^^) 




Wovelength {f^n^ 



Fig. 9. — Comparison between the best fitting solar-abundance 
model and the broad-band photometry (top panel) and OSIRIS 
broad- band H and K spectra (middle and bottom panels). Obser- 
vations are show in black with 1-cr error bars, the model spectrum 
in red and model photometry in green. The photometric upper 
limit from Hinz et al. (2010) is plotted as a horizontal bar. 



fit by minimization that samples a continuous distri- 
bution of parameters is beyond the scope of this paper. 
Instead, a was computed for each model spectrum 
from the grid compared to each of the three data sets 
(photometry, H and K-hand spectra). When computing 
a x'^ for the photometric comparison, each model spec- 
trum was convolved with the appropriate filter response 
functions to produce synthetic photometry in J, H, Kg, 
and L' bands along with the filters used by Hinz et al. 
(2010) and Currie et al. (2011); Mbar, z, and a filter cen- 
tered at 3.3 fim. When comparing models to the OSIRIS 
spectra, the synthetic spectra were first smoothed (with 
a Gaussian kernel) to the native instrumental response 
for OSIRIS (i? 4000), then binned down to the same 
wavelength sampling (and in the same manner) as done 
for the OSIRIS data cubes. The best fitting model was 
taken to be the model with smallest mean for the three 
data sets. Figure 8 shows the mean reduced-x^ distribu- 
tion across Tefi and log{g) for the model comparisons to 
all three sets of data. The best fit is located at Teff= 
HOOK and log(g)= 3.5. Formal error-bars for the fit 
have not been determined; however, models within the 
mean x^=4 contour are all reasonable fits to the data. 
Therefore, the error on the fit is close to ±100K for effec- 
tive temperature and ±0.5 dec for surface gravity, and 
comparable to most brown dwarf model fit uncertainties 
found in the literature. 

Figure 9 compares the best fit to each of the three 
data sets. The model spectrum agrees very well with the 
photometry. The differences seen at J, H, Kg, and L' 
are within the 1-cr error bars. The model only slightly 
under predicts the fiux at 3.3 /im, compared to the Currie 
et al. (2011) detection. The fiux measurement at 3.3 iim 
excludes a large number of models, especially at high 
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Toff, and is quite sensitive to the presence or absence of 
CII4. Observations of other systems (e.g., 2M1207b) in 
this band-pass would be very useful. 

At H band, the model spectrum is marginally consis- 
tent with the observed spectrum. The slope at the red 
and blue ends are not as steep as observed and, overall, 
the model spectrum is too fiat. It is possible that the 
residual flux from the bright speckle feature seen near ~ 
1.7/j,m (see Fig. 3) remains in the final observed spectrum 
and is responsible for the larger disagreement seen at this 
wavelength. As for the red and blue iJ-band slopes, these 
wavelength regions have historically been difficult to fit 
mainly due to missing CH4 line opacity. Consequently, 
the level of disagreement at i?-band is not too surpris- 
ing. At K band, the model compares favorably to the 
observed spectrum. The red and blue slopes are in good 
agreement as is the location of the peak flux. 

When the data sets are compared to the models indi- 
vidually, slightly different results are obtained. The best 
fitting model to the photometry has T^s = lOOOK and 
\og{g) = 5.0. H band is best fit by a model with T^s 
= 1200K and log{g) — 3.0 (clearly favoring low gravity); 
however, this model has colors that are too blue when 
compared to observed photometry, especially in the z- 
band. At K band, the best fit has T^s = HOOK and 
\og{g) = 3.0. When the mean x^ for H and K is min- 
imized, the best fit is TcS — 1200K and log(g) — 3.5. 
Each of these best fit models has a slightly different cloud 
thickness, but all have K^z — 10^ cm^ s~^ and ao = 5 
/im. 

5.2. Inferred Atmospheric Properties 

Despite the imperfect fit to the if-band, there is strong 
evidence for low surface gravity. In Figure 10 the i?-band 
OSIRIS spectrum is compared to a sequence of model 
spectra with different gravities but with all other param- 
eters equal to those of the overall best fit shown in Fig. 
9. The observed triangular shaped pseudo-continuum 
is best reproduced by the lowest-gravity model. The 
shape of the H band is often used as a gravity indica- 
tor in L-type brown dwarfs, where low gravity objects 
have distinct triangular spectra (Lucas et al. 2001; Kirk- 
patrick et al. 2006) compared to the more rounded spec- 
tra of older, higher gravity, field brown dwarfs. This 
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Fig. 11. — Same, best fitting model as shown in Fig. 9 (red), 
compared to models with the same parameters and temperature- 
pressure profile, but with CH4 opacities removed (blue) or both 
CO and CH4 opacities removed (green). 



gravity- dependence of the H-hand is produced by weak- 
ening of the continuous opacity produced by colhsion- 
induced absorption (CIA) as gravity decreases, relative 
to the molecular line opacity (most significantly relative 
to that of H2O). When CIA is weak, the iJ-band shape is 
defined by the natural triangular shape of the H2O bands 
on either side of ~ 1.7/im. However, when the surface 
gravity is high, the larger pressures increase the CIA to 
a point where it can become the dominate opacity source 
between the two H2O bands on either side of the peak. 
The smooth wavelength dependence of CIA rounds the 
i?-band peak (Borysow et al. 1997; Kirkpatrick et al. 
2006; Rice et al. 2011). For field brown dwarfs, the effec- 
tiveness of iJ-band as gravity indicator is decreased for 
field T dwarfs as CH4 opacity dominates the near-IR and 
competes with that of II2O around the H band. How- 
ever, if CH4 mole fractions are reduced by vertical mix- 
ing and non-equilibrium chemistry at low gravity and low 
Teff, as is likely the case for HR8799b, then the behavior 
of the H band with gravity seen in hotter L-dwarfs will 
persist even down to effective temperatures most often 
associated with T-dwarfs. 

Water is by far the main molecular opacity source shap- 
ing the OSIRIS spectra. The combination of high K^z 
and low \og{g) leads to, in the models used here, a sharp 
reduction in the CH4 mixing ratio (e.g. see Fig. 7) with 
a correspondingly large increase in the CO mixing ra- 
tio. Figure 11 compares the best fit model to those with 
the same parameters and temperature-pressure struc- 
ture, but with CO or CH4 opacities removed. Weak CH4 
absorption is evident in the K band, as is a possible con- 
tribution from CO; however, the strongest evidence for 
CH4 absorption is the Currie et al. (2011) detection (and 
Hinz et al. (2010) upper limit) at 3.3 /im where the flux is 
clearly suppressed by CH4. Unfortunately, the Hinz et al. 
(2010) upper limit near 4.5 fim is not sensitive enough to 



distinguish between the majority of the models explored 
here and, therefore, does not provide strong constraints 
on CO. 

The best fitting, solar abundance, atmosphere model 
for HR8799b (shown in Fig. 9) has a cloud base just be- 
low P = 1 bar with a composition that is typical for 
brown dwarf models with Teff > lOOOK. Magnesium- 
silicate (e.g. forsterite) and iron grains are the most 
abundant, with a mixture of other grains with lower mole 
fractions. The structure of the cloud has many charac- 
teristics in common with brown dwarf cloud predictions, 
namely an abrupt drop in the gas-to-dust ratio at higher 
pressures and temperatures, where condensation is no 
longer thermodynamically favored. This cloud base is 
determined by the intersection of the T-P profile and 
the relevant condensation curves (see Fig. 12). At lower 
temperatures and pressures, the gas-to-dust ratio drops 
exponentially (by design in this study, see lower panel 
of Fig. 12) in a similar manner as the smooth decrease 
in dust predicted by the models of Helling & Woitke and 
Marley & Ackerman (see Helling et al. (2008a) for a com- 
parison). Location and vertical extent are the most im- 
portant aspects of the cloud structure for HR8799b pre- 
ferred by this work. Generally, the cloud base retreats 
to larger pressures with decreasing Toff and fixed grav- 
ity; however, the low gravity of the best fitting HR8799b 
model counteracts this trend. As gravity decreases the 
T-P profile shifts to lower pressure, moving the cloud 
base along with it. At logig) ^ 3.5, the cloud base sits 
near the photosphere with vertical extent that is far less 
than an equilibrium cloud model (see Fig. 12), but suf- 
ficient to produce redder colors than cloud-free brown 
dwarf models with the same Tes and \og{g). 

The modal particle size for HR8799b, favored by this 
work is around 5 to 10 fim. Grains of this size are 
within the range predicted by the Marley & Ackerman 
and Helling & Woitke for models with Teff - lOOOK. 
Both groups predict sizes between 1 and 100 /zm at high 
pressures (> 10 bar) and sizes that decrease below a mi- 
cron at low pressures (< 1 bar). Burrows et al. (2006) 
chose a.0— 100 /im (at the high end of most predictions) 
for their set of reference model atmospheres arguing that, 
in general, L-type brown dwarfs are best fit when larger 
grain sizes are assumed but acknowledged that such large 
grains make fitting the reddest late L dwarfs difficult. As 
already discussed above, the observed near-IR spectra of 
the reddest late L and early T dwarfs come closest to 
matching the H and K-hand spectra of HR8799b and, 
thus, perhaps it should not be surprising that the best 
model fits are found among those with smaller grains. 
However, the simulations of Cooper et al. (2003) predict 
that particle size should actually increase with decreasing 
gravity. When comparing grain sizes adopted by differ- 
ent models, it is important to remember that grain size 
alone does not determine the total cloud opacity - com- 
position also plays an important role. In the cloud model 
used here a mixture of grains are included in the total 
opacity and such a cloud can be more opaque than a 
cloud including opacity from only pure silicate grains of 
the same size. 
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Fig. 12. — Top : Temperature versus pressure for the best fit- 
ting model identified in Fig. 8 (solid line) and for two models with 
the same T(,ff, log((ji) and Kzz, but with different cloud model pa- 
rameters (dotted lines). Also shown (dashed line) is the best fitting 
metal rich model (Fig. 14). Condensation curves (thick solid lines) 
for Fe and MgSiOs are plotted to indicate the general location of 
the cloud base, located at the intersection of these curves. Other 
grains (e.g. Mg2Si04, AI2O3, Ni, etc.) have condensation curves 
in same general location. Bottom : Dust fraction versus pressure 
for the same models shown in the top panel (in order from left to 
right). A constant dust fraction with decreasing pressure is a hall- 
mark of the pure chemical equilibrium cloud (e.g. the 'DUSTY' 
models of AUard et al. 2001). In contrast, the intermediate cloud 
model (ICM) used here can have lower peak dust fraction than the 
equilibrium cloud and decays with height (see text for details). 



Since the distance to HR8799 has been precisely mea- 
sured (39.4 ± 1.0 pc, van Leeuwen 2007), the luminosity 
of HR8799b is known (logLobs/-^© = —5.1) and, there- 
fore, radii can be assigned to each Teg value on the left 
vertical axis of Fig. 8; these radii are given by the right, 
non-linear, vertical axis. Using these radii and the log (17) 
axis values, lines of constant mass can be plotted that are 
completely independent of interior cooling track models. 
Dynamical stability arguments suggest that, in order for 
the three planets to be as massive as indicated by the 
cooling tracks, they must be in double 2:1 resonances 
and have masses < 20 Mjup (Fabrycky & Murray-Clay 
2010). Using 20Mjup as an upper limit on the mass of 
HR8799b sets a Toff/log(g) boundary in Fig. 8 (right- 
most solid line). This limit was the motivation for not 
extending the model grid to log(f/) greater than 5.0. A 
second boundary, at lower surface gravity, can be set 
using the Kelvin- Helmholtz time scale (tkh)i which re- 
flects the approximate scaling of detailed planetary cool- 
ing curves. Since the planets almost certainly are older 
than 1 Myr, tkh = 1 Myr is a very conservative limit 
and is indicated by the left-most dotted line in Fig. 8. 
In other words, planets located to the left of this dot- 
ted line could not maintain the observed luminosity of 
HR8799b for more than 1 Myr (under the usual assump- 
tion that the primary energy source is gravitational po- 
tential energy). This time scale limit was used to exclude 
a second minimum at log(5) = 2.5 near Tcff = 1100 and 
1200K. For Teg = HOOK ±100K, the radius must be 0.75 
R'Jupio'12 to match Lobs- Using this radius and the best 
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Fig. 13. — Plotted are lines of constant mass (solid) and 
isochrones (dotted) from Baraffe et al. (2003) cond evolution mod- 
els. The shaded region covers the range of T^fr and log(g) predicted 
by these evolution models for an object with the same luminosity 
as HR8799b (the darker region brackets 10 and 120 Myrs). The 
best fitting Teg and log((;), based on fitting solar abundance atmo- 
sphere models (see text for details) , is indicted by the filled symbol 
surrounded by the = 4 contour from Fig. 8. 



fit surface gravity, \o'g{g) = 3.5 ± 0.5, the corresponding 
mass is just under 1 Jupiter-mass (~ 0.72 Mjuplo;g). a 
factor of ~ 10 lower than predicted by evolutionary cool- 
ing tracks. A radius of 0.75 Rjup is also ^ 30% smaller 
than cooling track predictions. 

While the implied mass of ^ 0.7 Mjup is comfortably 
within the ranged needed for dynamical stability, it is dif- 
ficult to form objects, either by core- accretion or by grav- 
itational instability, with this mass, radius (~ 0.75 Rjup), 
and the planet's observed luminosity. This combination 
is inconsistent with our understanding of degenerate or 
partially-degenerate H+Hc gaseous bodies with any rea- 
sonable heavy-metal composition or of any age (Baraffe 
et al. 2008). Figure 13 compares the Toff and log(g) val- 
ues from the atmosphere model fit to those predicted by 
evolution models from Baraffe et al. (2003). While the 
best fit surface gravity is within the expected range, an 
effective temperature of HOOK is ~ 300K above the tem- 
peratures predicted by all current hot-star cooling tracks, 
for an age less than 100 Myrs. The high effective tem- 
perature would be consistent with evolution predictions 
if the system was old (> 1 Gyr) and log((7) was greater 
than 5. However, as shown above, such a high surface 
gravity would be very difficult to reconcile with the H- 
band spectrum and would have serious dynamical stabil- 
ity implications (Fabrycky & Murray-Clay 2010; Marois 
et al. 2010; Currie et al. 2011). 

Giant planet evolution models are known to become 
increasingly unreliable at young ages for the simple rea- 
son that planet formation is not well understood and, 
thus, the initial conditions are not known. This fact was 
highlighted by Baraffe et al. (2002), with a focus on the 
formation of brown dwarfs and low mass stars from the 
collapse of molecular clouds (the main realm of appli- 
cability for the common cond and dusty cooling tracks). 
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Fig. 14. — Same, best fitting, model as sliown in Fig. 9 (red), 
compared to tlie best fitting model that lias T^g and log{g) con- 
sistent with standard "hot start" evolution models (green). Also 
plotted are the synthetic spectra and photometry for the chemical 
equilibrium atmosphere model (blue) from Marois et al. (2008). 
See text for details. 



More recently, Marley et al. (2007) revisited the question 
of initial conditions by comparing the evolution of plan- 
ets starting with high Tes and luminosity ( "hot start" , 
as often adopted for brown dwarfs and motivated by the 
Stabler birth line), to the evolution of planets starting 
with the colder and less luminous initial conditions ( "cold 
start" ) that emerge from core-accretion planet formation 
models. Marley et al. concluded that young planets 
that formed by core-accretion are cooler, smaller, and 
thus fainter than planets formed by single-mode gravi- 
tational collapse within a disk or giant molecular cloud. 
Given that HR8799b is clearly underluminous compared 
to the hot-start models and overluminous for the cold- 
start models (Marley et al. 2007; Fortney et al. 2008), 
it is unlikely that initial conditions alone are responsible 
for the offset between the atmosphere-inferred properties 
and the evolution track predictions (a situation similar 
to that of 2M1207b, as discussed below). 

A ~ 300K offset in effective temperature is similar in 
magnitude to offsets in effective temperature found for 
a number of M, L and T-type binary brown dwarfs (Liu 
et al. 2008, 2010; Konopacky et al. 2010; Dupuy et al. 
2010). This suggests that a systematic offset may exist 
between the atmosphere and interior model predictions. 
A 10 to 20 % difference in TcS is comparable to the dif- 
ferences between cloudy and cloud-free tracks in the 10 
and 100 Myr range and 10 Mjup(Chabrier et al. 2000), 
so it is still unclear what fraction of such an offset could 
be resolved by improvements purely in the modeling of 
the interior or the atmosphere. However, the increas- 
ing number of free parameters in substellar atmosphere 
models is a clear indication that they are even more com- 
plex than previously imagined. Thus, the possibility still 
remains that hidden among the free parameters is an at- 



mosphere model with Toff and log(g) that matches the 
evolution model predictions. 

6.1. Non-solar Abundance Model Comparisons 

The need for an intermediate cloud model (constrained 
vertically), to match the observed broad-band photom- 
etry of all three HR8799 planets, was first discussed by 
Marois et al. (2008). Marois et al. (2008) showed that 
a synthetic spectrum from a solar abundance model at- 
mosphere with an intermediately thick cloud and Toff / 
\og{g) consistent with evolution model predictions (820K 
/ 4.0) matched the available observed photometry at 3- 
(7. While only CH4s//, J, H, Kg, and L' photometry 
where available in 2008, this earlier (chemical equilib- 
rium) model remains consistent (^ S-cr) with all of the 
revised/new broad-band photometry (see Fig. 14). How- 
ever, the true level of agreement between this model and 
reality is revealed by the OSIRIS data (Fig. 14) where 
much weaker CH4 absorption is observed than predicted 
by this chemical equilibrium model. This discrepancy 
highlights the difficulties in obtaining a reliable descrip- 
tion of young giant planets from photometry alone, which 
can hide significant clues about their atmospheric prop- 
erties. 

To explore differences between the atmosphere and 
evolution derived bulk parameter more closely, a sepa- 
rate grid of atmosphere models was computed, with Tcff 
and log{g) uniformly sampling the evolution track pre- 
dictions that are consistent with the observed luminos- 
ity of HR8799b (down the center of the grey band in 
Fig. 13). The clouds and non-equilibrium chemistry were 
modeled in the same manner as in the solar abundance 
grid described above, yet only ao = 5 and 10 fim were 
used. The metallicity, however, was extended beyond so- 
lar abundances to include two metal rich cases ([Fe/H] 
= 0.5 and 1.0) and an abundance pattern matching the 
observed peculiar (A-boo-type) subsolar abundances of 
HR8799 (Sadakane 2006). The best fitting model from 
this "evolution-consistent" grid has Tcff = 896K, log(g) 
= 4.3, [Fe/H] = 1.0, logXboi/i© = -5.06 (see Fig. 14). 
The agreement between this metal-rich model and the 
OSIRIS H and /C-band spectra is similar to that of the 
previous hotter, solar abundance, best fit. The broad- 
band photometric comparison is also good (mostly to 1- 
fj), but slightly worse compared to the solar abundance 
model at J and K bands (see Figure 14) but still consis- 
tent at 2-a. The lowering of Tgft by increasing metallic- 
ity is to be expected. Looking at the Mj versus J — K 
CMD in Fig. 6, in order for cooler atmospheres to match 
HR8799b they need to be redder in J — and this can 
be accomplished by increasing metallicity and/or by de- 
creasing the modal grain size (Burrows et al. 2006). The 
mean for this metal-rich comparison is ~ 5.5 and only 
marginally better than the solar abundance model at the 
same temperature and gravity. However, given that even 
in this smaller grid the parameter space is still sparsely 
sampled, there is plenty of room for improvement. 

6.2. Comparison to 2M1207 b 

In near-IR CMDs, the object most similar to HR8799b 
is 2M1207b, a planet-mass companion to a young 10 
Myr) brown dwarf (Chauvin et al. 2005). In addition 
to their very red colors, these two objects also share the 
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Fig. 15.— OSIRIS H and X-band spectra of HR8799b compared 
to the spectra of 2M1207b (Patience et al. 2010) binned to the same 
resolution and interpolated onto the same wavelength grid as the 
OSIRIS observations. The 2M1207b spectra where independently 
normalized for this comparison. The error-bars for the 2M1207b 
spectra are significantly smaller than those of HR8799b and, for 
clarify, are not show. 



same puzzling fact that their SEDs appear hotter than 
evolution models predict, for their ages and luminosi- 
ties. Available photometry and spectroscopy of 2M1207b 
show all the signs of an object with effective temperature 
~ 1600K (Mohanty et al. 2007; Ducourant et al. 2008; 
Patience et al. 2010). When combined with current lumi- 
nosity estimates (log(L/L0) ~ —4.5), the radius implied 
by this high effective temperature is 0.7 Rjup, com- 
parable to the small size inferred above for HR8799b. 
Consequently, both 2M1207b and HR8799b are under- 
luminous compared to hot-start evolution models and 
overluminous for cold start models, for their respective 
ages. 

Mohanty et al. (2007) concluded that the high effective 
temperature was likely correct and, thus, 2M1207b is un- 
derluminous by ~ 2.5 magnitudes. Given that the dis- 
tance to the source (52.4 ±1.1 pc) is known (Ducourant 
et al. 2008), Mohanty et al. (2007) proposed that an edge- 
on disk surrounds 2M1207b providing the necessary ex- 
tinction. Alternatively, Mamajek & Meyer (2007) pro- 
posed that 2M1207b is not a normal gas giant but rather 
the product of a recent protoplanetary collision and the 
observed flux is the afterglow from the hot, physically 
small, remnant. Given that HR8799b suffers from a sim- 
ilar quandary, it is unlikely that such rare events should 
be necessary to explain the first and only two objects 
found in the very red and faint regions of near-IR CMDs. 

Given the puzzling similarities between 2M1207b and 
HR8799b mentioned above, it is useful to compare their 
SEDs. Moderate resolution {R - 500 - 1500) H and 
if-band spectra of 2M1207b (Patience et al. 2010) are 
compared to the HR8799b spectra in Figure 15. The 
spectrum of 2M1207b has shallower H2O bands with a 
distinct CO band-head at 2.3 fim and no prominent CH4 
absorption - a spectrum very much like typical L dwarfs. 
The only feature the two sets of spectra have in com- 
mon is a triangular shaped spectrum in the H band, 
consistent with both objects having low surface gravity 
(indicative of low mass and/or youth). Assuming both 
objects are not the result of recent protoplanetary colli- 
sions, 2M1207b and HR8799b may well represent impor- 
tant evolutionary states of substellar atmospheres. In- 



dependent of the edge-on disk hypothesis, if the ages of 
the two systems are correct, then the planets' masses are 
probably the same to within a factor of 2. If the masses 
are similar, then these objects provide evidence for rapid 
spectral changes from 10 to ^ 30 Myr for young planet- 
mass objects. 

The analysis of HR8799b described above provides a 
third possible explanation for the properties of 2M1207b 
- namely that certain combinations of cloud coverage 
and non-equilibrium chemistry (and perhaps metallic- 
ity) are capable of producing a very low temperature 
object with an L-type near-IR spectrum. The two hy- 
potheses proposed by Mohanty et al. and Mamajek et 
al. hinge on the assumption that the effective temper- 
ature of 2M1207b is ^ 1600K. However, looking back 
at Figure 6, 2M1207b is intersected by an intermedi- 
ate cloudy model with T^g = HOOK, close to evolution 
model predictions and consistent with a typical radius. 
Also, given the apparently extreme reduction of CH4 in 
the atmosphere of HR8799b, it is conceivable that the 
CO/CH4 mole fractions are quenched at even higher 
pressures in the younger (and likely lower gravity) at- 
mosphere of 2M1207b, resulting in a low temperature 
atmosphere that is rich in CO (perhaps enough to pro- 
duce a distinct CO bandhead at 2.3 /im, as seen in the 
spectrum from Patience et al.) and without significant 
amounts of CH4. A search for such a cold, cloudy, model 
atmosphere with extreme non-equilibrium chemistry for 
2M1207b is underway and will be presented in a separate 
paper (Barman et al., in prep.). 

Despite their spectral differences, the most straight- 
forward explanation for the positions of 2M1207b and 
IIR8799b in near-IR CMDs is one of atmospheric ori- 
gins. Instead of an edge-on disk or recent protoplanetary 
collision, a full exploration of the cloud properties, non- 
equilibrium chemistry, and metallicity is likely required 
to reproduce their observed SEDs with a model atmo- 
sphere having T^g and \og{g) consistent with evolution 
predictions. However, as more and more important phys- 
ical processes are identified in substellar atmospheres, 
the number of free (and essentially independent) param- 
eters increases and the ability to exhaustively explore all 
possible parameter combinations diminishes - more work 
(along the lines of Freytag et al. 2010) is clearly required 
to reduce the current set of parameters and produce more 
unified model atmospheres. 

7. CONCLUSIONS 

The HR8799 planetary system provides a rare oppor- 
tunity to study multiple, coeval, giant planets at very 
young ages. The data and model comparisons described 
above are only first steps toward understanding just one 
of these planets and many challenging problems remain. 
While the data presented here offer the broadest spectral 
coverage for an HR8799 planet obtained so far, the low 
resolution and modest signal-to-noise set strong limits on 
what can be understood about HR8799b. In addition to 
the limits imposed by the data, the complexity and excess 
of free parameters in model atmospheres for giant planets 
and brown dwarfs also set strong limits on what can be 
confidently inferred from observations. With these lim- 
itations in mind, several conclusions can be drawn from 
the observations of HR8799b. 

The SED of HR8799b is shaped by water absorp- 
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tion bands that arc deeper than typieally seen in late 
L and early T field dwarfs, but clearly eonsistent with 
a Hydrogen-rich planetary atmosphere. Only weak 
methane absorption is detected in the K band along 
with very weak CO absorption, indicating strong non- 
equilibrium chemistry and vertical mixing. Additional 
evidence for very weak methane absorption is seen in 
the H band spectrum (and narrow-band photometry). 
The near-IR spectra and broad-band photometry are 
smoothed and reddened by photospheric clouds not typ- 
ical for brown dwarfs with Teff< HOOK. The triangular 
shape of the iJ-band spectrum clearly indicates low sur- 
face gravity (log(g) < 5), consistent with an age signifi- 
cantly yoimger than that proposed by Moya et al. (2010). 

The best match from a solar abundance grid of model 
atmospheres (having Toff= HOOK and log(g)= 3.5) 
found above implies a radius that is too small to be con- 
sistent with both the observed luminosity and the current 
understanding of the interiors and evolution of H-|-He- 
rich planets, regardless of mass, age, or interior core-size. 
Furthermore, a variety of model atmospheres can be con- 
structed that fit the observed broad-band photometry 
very well (at least from the model grids used here), yet 
do not match the observed spectra - indicating that pa- 
rameters inferred from broad-band photometric fits alone 
should be treated with caution. 

While certainly not perfect, the model comparisons 
presented here are better representations of the obser- 
vations than those presented by Bowler et al. (2010) 
and Currie et al. (2011). In particular. Bowler et al. 
found fits to their narrow-band spectrum that require 
Tefi> HOOK, which leads to unphysically small radii of 
R < 0.75Rjup. These previous studies also had diffi- 
culty simultaneously fitting the photometry; for exam- 
ple, Currie et al. found it difficult to match most of 
the photometry, especially at A > 3/Ltm (even with their 
thick or patchy cloud models) and Bowler et al. had dif- 
ficulty matching the near-IR colors for reasonable Tcff. 
These difficulties may stem from the sparse sampling of 
model parameters in the grids used in these two stud- 
ies (e.g. those parameters regulating non-eq. chem- 
istry and clouds), though considerable degeneracies are 
present when fitting only broad-band photometry. Fur- 
thermore, the best fitting models for HR8799b found by 
Currie et al. typically have deeper methane absorption 
than seen in the broad K-h&nA spectrum described above 
or the narrow-band spectrum of Bowler et al. Interest- 
ingly, Currie et al. emphasize cloud thickness/coverage 
over non-equilibrium chemistry, in particular as it per- 
tains to the near-IR SED. The analysis presented here, 
however, requires deep quenching of CO/CH4 to match 
the near-IR spectra (and near-IR narrow-band photom- 
etry), placing equal importance on both atmospheric ef- 
fects when reproducing the observations. 

Finding a solar abundance atmosphere that matches 
the complete set of data and interior /cooling-track pre- 
dictions seems unlikely (though still not inconceivable). 
A more likely scenario is that the atmosphere of IIR8799b 
has an elemental composition different from the Sun. A 
model atmosphere that matches the observations reason- 
ably well (as well or better than any previous study) , in- 
cluding the predictions of evolution models, can be found 
if the atmosphere is enhanced by a factor 10 in metals 
compared to the Sun. Such a model has Teg = 896K and 



\og{g) = 4.3, consistent with a young giant planet with 
mass less than 10 Mjup. At ^ 30 Myrs, as recently pro- 
posed by Zuckerman et al. (2011), the mass of HR8799b 
would be 4 to 5 Mjup. 

If HR8799b is indeed metal rich and near the very low 
end of the original mass range inferred by Marois et al. 
(2008), then core-accretion may be the preferred forma- 
tion mechanism for the HR8799 planetary system. Metal 
enrichment, relative to the host star, is often considered 
a hallmark of the core-accretion scenario and recent work 
indicates that if the HR8799 planets formed by gravita- 
tion instabilities, they should show little metal enhance- 
ment (Helled & Bodenheimer 2010). Atmosphere models 
with either solar abundances or the metal poor A-boo-like 
abundances of the star result in Tefi that are too high and 
corresponding radii that are too small. The analysis of 
Bowler et al. (2010) also indicates that enhanced metals 
(and low gravity) are likely to produce the best model 
fits. 

While these conclusions rest heavily on models not 
tested by a large population of young giant planets 
(such a population should emerge in the coming years), 
the conclusions are self-supporting. Low surface grav- 
ity, clearly preferred by the model fitting (regardless 
of metallicity) and indicated by the shape of the H- 
band spectrum, may well allow for more efficient ver- 
tical mixing. As surface gravity decreases, the radiative- 
convection boundary moves closer to the photosphere, 
the maximum convective velocity increases, and K^z also 
increases, at least in the convection zone (Freytag et al. 
2010; Rice et al. 2011). Enhanced vertical mixing leads 
to enhanced non-equilibrium chemistry, as suggested by 
the absence of CH4 absorption in both H and K bands. 
However, the connection between gravity and mixing is 
merely suggestive as the primary physical mechanism for 
mixing into the radiative zone unknown. Low surface 
gravity may also establish atmospheric conditions (not 
found in the higher gravity field T dwarfs) that support 
thick clouds at photospheric depths, allowing the inter- 
section of the T-P curve and the condensation curve to 
be near the photosphere even at low Teff. A connection 
between clouds and surface gravity is already starting 
to emerge for brown dwarfs. Evidence for a gravity de- 
pendence for the L-to-T transition has already been re- 
ported (Metchev & Hillenbrand 2006; Saumon & Marley 
2008) and, recently, Stephens et al. (2009) found that 
low-gravity brown dwarfs in their sample remain cloudy 
at lower effective temperatures, with the L-to-T tran- 
sition temperature decreasing from 1300K to HOOK as 
log{g) decreases from 5.0 to 4.5. Metchev & Hillenbrand 
(2006) also found that HD20303B, a low gravity tran- 
sition brown dwarf, is ^ 250K cooler than field objects 
with similar spectral types. If HR8799b has a log(g) 
~ 4, then it is possible that this emerging trend ex- 
tends to this object as well, and would suggest that very 
low gravity objects (like young giant planets) are cloudy 
down to temperatures well below lOOOK. The structure 
of the clouds needed to reproduce the colors and spectra 
of the HR8799 planets and 2ml207b are similar to those 
needed for late L-dwarfs. This suggests that such ob- 
jects, from an atmospheric point of view, are extensions 
of the L-type field dwarfs rather than members of a new 
class altogether. This conclusion is distinct from that of 
Currie et al. who claim the clouds must be significantly 
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thicker than those required to match L dwarfs. Finally, 
enhanced metals would also lead to enhanced grain for- 
mation and conditions ripe for atmospheric cloud forma- 
tion. The combination of these atmospheric conditions 
lead to bulk properties, inferred from the SED analysis 
alone, that include a radius consistent with the observed 
luminosity and current understanding of giant planet in- 
teriors, as well as a mass that fits comfortably within the 
limits required for long-term dynamical stability. 
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